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ABSTRACT 

We report on a multiwavelength observation of the blazar 3C 454.3 (which 
we dubbed crazy diam,ond) carried out on November 2007 by means of the astro- 
physical satelhtes AGILE, INTEGRAL, Swift, the WEBT Consortium, and the 
optical-NIR telescope REM. Thanks to the wide field of view (FoV) of the AGILE 
satellite and its prompt alert dissemination to other observatories, we obtained a 
long (three weeks), almost continuous 7-ray coverage of the blazar 3C 454.3 across 
fourteen decades of energy. This broad-band monitoring allows us to study in 
great detail light curves, correlations, time-lags and spectral energy distributions 
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(SEDs) during different pliysical states. Gamma-ray data were collected during 
an AGILE pointing towards tlie Cygnus Region. Target of Opportunity (ToO) 
observations were performed to follow-up the 7-ray observations in the soft and 
hard X-ray energy bands. Optical data were acquired continuously by means of 
a pre-planned WEBT campaign and through a REM ToO repointing. 3C 454.3 
is detected at a ~ 19 — a level during the 3-week observing period, with an aver- 
age flux above 100 MeV of FE>iooMeV = (170 ± 13) x 10"^ photons cm'^ s'^. 
The 7-ray spectrum can be fit with a single power-law with photon index 
Egrid = 1-73 ± 0.16 between 100 MeV and 1 GeV. We detect significant day- 
by-day variability of the 7-ray emission during our observations, and we can 
exclude that the fluxes are constant at the 99.6% (~ 2.9a) level. The source 
was detected typically around 40 degrees off-axis, and it was substantially off- 
axis in the field of view of the AGILE hard X-ray imager. However, a 5-day 
long ToO observation by INTEGRAL detected 3C 454.3 at an average flux of 
about F2o-2ookeV = 1-49 X 10"'^ photons cm~^ s~^with an average photon index of 
FiBis = 1.75 ± 0.24 between 20-200 keV. Swift also detected 3C 454.3 with a flux 
in the 0.3-10 keV energy band in the range (1.23 — 1.40) x 10^^ photons cm^^ s~^ 
and a photon index in the range Fxrt = 1-56 — 1.73. In the optical band, both 
WEBT and REM show an extremely variable behavior in the R band. A correla- 
tion analysis based on the entire data set is consistent with no time-lags between 
the 7-ray and the optical flux variations. Our simultaneous multifrequency obser- 
vations strongly indicate that the dominant emission mechanism between 30 MeV 
and 30 GeV is dominated by inverse Compton scattering of relativistic electrons 
in the jet on the external photons from the broad line region. 

Subject headings: galaxies: active - galaxies: quasars: general - galaxies: 
quasars: individual: 3C 454.3 - galaxies: jets - radiation mechanism: non ther- 
mal 



Introduction 



Among activ e galactic nucle i (AG Ns), blazars show intense and variable 7-ray emission 
above 100 MeV (IHartman et al.lll999l ). Variability timescale can be as short as few days, 
or last few weeks. They emit across several decades of energy, from the radio to the TeV 
energy band. Blazar spectral energy distributions (SEDs) are typically double humped with 
a flrst peak occurring in the IR/Optical band in the so-called red blazars (including Flat 
Spectrum Radio Quasars, FSRQs, and Low-energy peaked BL Lacs, LBLs) and at UV/X- 
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rays in the so-called blue hlazars (including High-energy peaked BL Lacs, HBLs) and it is 
commonly interpreted as synchrotron radiation from high-energy electrons in a relativistic 
jet. The second SED component, which peaks at MeV-GeV energies in red blazars and at 
TeV energies in blue blazars, is commonly interpreted as Inverse Compton (IC) scattering of 
soft seed photons by relati vistic electrons. A re c ent re view of the blazar emission mechanisms 
and energetics is given in lCelotti &: Ghisellinil (120081 ). Alternatively, the blazar SED can be 
explained in the framework of the hadronic models, where the relativitic protons in the jet 
are the primary accelerate d particles, emitting 7-ray radiation by mean s of photo-pair and 
photo-pion production (see lMiicke &: Protheroell200ll ; iMiicke et al.ll2003l for a recent review 
on hadronic models). 

Multiwavelength studies of variable 7-ray blazars have been carried out since the be- 
ginning of the 1990s, thanks to the CGRO Observatory. Nevertheless, only a few objects 
were detected on a timescale of about two weeks in the 7-ray energy band, and simulta- 
neously monitored at different energies, obtaining a multifrequency coverage. Among the 
FSRQs detected at energies above 100 MeV by th e EGRET telescope on-board the Compton 
Gamma-Ray Observatory jHartman et^ll999l ). 3C 454.3 (PKS 2251+158; z = 0.859) is 
certainly one of the most active at high energy. In the EGRET era , it was detected in 1992 
during an intense 7-ray flaring episode (jHartman et al.lll992l . Il993l ) when its flux -FE>iooMeV 



was observed to vary within the range (0.4 — 1.4) x 10 ^ photons c m ^s ^. In 1995, a 2-week 



campaign detected a 7-ray flux < 1/5 of its historical maximum (lAUer et al.l 119971 ). 



In 2005, 3C 45 4.3 underwent a major flaring activity in alm ost all energy bands (see 
Giommi et al.l l2006l ). In the optical, it reached R = 12.0 mag (jVillata et al.l l2006l ) and 



it was detected by I NTEGRAL at a level of ~ 3 X 10 ^photons cm ^ s ^ in the 3- 



200 keV energy band (jPian et al.ll2006l ). Since the detection of the exceptional 2005 outburst, 
several monitoring campaigns were carried out to fo llow the source multifrequency behavior 
(IVillata et aPbood . I2OO7I : iRaiteri et allboO?! . l2008al lbll During the last of these campaigns, 
3C 454.3 underwent a new optical brightening in mid July 2007, which triggered observations 
at all frequencies. 



In July 2007, IVercellone et al.l (l2008l . hereafter V08) reported the highest 7-ray flare 
from 3C 454.3. During the period 2007 July 24-30 the average 7-ray flux was -FE>iooMeV = 
(280 ±40) X 10^® photons cm~^ s~^, about a factor of two higher than in 1995. No emission 
was detected by Super-AGILE in the energy range 20-60 keV, with a 3-0" upper limit of 
2.3 X 10~^ photons cm~^ s~^. 



By means of AGILE preliminary flux estimate (IVercellone et al.l 120071 ). iGhisellini et al. 



^Assuming a Crab- like spectrum. 
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( 120071 ) compared the 3C 454.3 SEDs obtained during three multiwavelength campaigns 
(2000, 2005 and 2007). The 2007 data show that the 7-ray high state occurred during a 
weaker optical/X-ray flux compared with the 2005 flare. 

In this paper (Paper I) we present the results of a multiwavelength campaign on 3C 454.3 
during a long-lasting 7-ray activity period between 20 7 Nov ember 10 and December 1. Pre- 
liminary 7-ray results were distributed inlChen et al.l (120071 ). while radio-to-optical and UV 
data were published in iRaiteri et al.l (j2008al ). A companion paper (Paper II, Donnarumma 
et al., in preparation) will describe the AGILE multiwavelength campaign during December 
2007. In Sect. [2] we present the multiwavelength campaign on 3C 454.3; in Sect. [3] through [6] 
we present the AGILE/GRID, INTEGRAL/IBIS, Swzft/XKT, WEBT and REM data anal- 
ysis, respectively; in Sec. [7] we present the simultaneous multiwavelength light-curves and 
SEDs. In Sect. [S] and [3 we discuss our results and draw our conclusions. Throughout this 
paper the quoted uncertainties are given at the l-cr level, unless otherwise stated. 



2. The Multiwavelength Campaign 



In 2007 November AGILE began pointing 3C 454.3 at high off-axis angle (about 40°) . 
Nevertheless, in a few days 3C 454.3 was dete cted at more than 5-cr ( jChen et al.l 120071 ). 
exhibiting variable activity on a day time-scale (jPucella et al.l 120071 ). Immediately after the 
source detection, a multiwavelength campaign started. AGILE data were collected during 
two different periods, the first ranging between 2007-11-10 12:17 UT and 2007-11-25 10:57 UT 
and the second between 2007-11-28 12:05 UT and 2007-12-01 11:39 UT, for a total of about 
592 ks. The three-day gap between them was due to a pre-planned GRID calibration activity. 
INTEGRAL data were collected during a dedicated ToO on revolutions 623 (between 2007- 
11-20 03:35 UT and 2007-11-22 08:46 UT) and 624 (between 2007-11-22 20:45 UT and 2007- 
11-24 15:50 UT), for a total of about 300 ks, while Swift/XKT data were obtained during 
several ToO pointings for a total of about 10 ks. WEBT data (radio to optical) as well as 
Swift/UVOT data were published in iRaiteri et al.l (l2008al ). while REM data were acquired 
following a ToO request. In both cases, optical data were acquired continuously during the 
whole AGILE campaign. 
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3. AGILE data 
3.1. Data Reduction and Analysis 



The AGILE satellite jTavani et aPboOSal lbl). a mission of the Italian Space Agency (ASl) 



devoted to high-energy astrophysics, is currently the only space mission capable of observing 
cosmic sources simultaneously in the energy bands 18-60 keV and 30 MeV-50 GeV. The 
satellite was launched on 2007 April 23 by the Indian PSLV-C8 rocket from the Satish 
Dhawan Space Center SHAR, Sriharikota. 

The AGILE scientific Instrument is very compact and combines four active detectors 
yielding broad-b and coverage from hard X-rays to gamma-rays: a Silicon Trac ker (ST; 



2007 



Prest et all2003l . 30 MeV-50 GeV), a co-aligned coded-mask hard X-ray imager (S A; 



100 MeV), and a segmented Anti-Coincidence System (ACS; Perotti et al. 20061 ). Gamma- 



Feroci et al. 



18-60 keV), a non-imaging Csl Mini-Calorimeter (M CALjabanti et allbood . 0.3 



ray detection is obtained by the combination of ST, MCAL and ACS; these three detectors 
form the AGILE Gamma- Ray Imaging Detector (GRID). 

Level-1 AGILE-GRID data were analyzed using the AGILE Standard Analysis Pipeline 
(see V08 for a detailed discussion of the AGILE data reduction). Since 3C 454.3 was at high 
off-axis angle, an ad-hoc 7-ray analysis was performed. We used 7-ray events filtered by 
means of the FT3ab_2 AGILE Filter pipeline. Counts, exposure, and Galactic background 
7-ray maps are created with a bin-size of 0?25 x 0?25 , for E > 100 MeV. Since the source 
was at 40° off-axis, all the maps were generated including all events collected up to 60° 
off-axis. We rejected all the 7-ray events whose reconstructed directions form angles with 
the satellite-Earth vector smaller than 80° (albrad=80), reducing the 7-ray Earth albedo 
contamination by excluding regions within ~ 10° from the Earth limb. The most recent 
versions (BUILD-15) of the Calibration files, which will be publicly availa ble at the ASI 
Science Data Centre (ASDC@) site, and of the 7-ray diffuse emission model (iGiuliani et al.l 



20041 ) were used. The first step consists in running the AGILE Source Location task in order 
to derive the most plausible location of the source. In a second step, we ran the AGILE 
Maximum Likelihood Analysis (ALIKE) using a radius of analysis of 10° , and the best guess 
position derived in the first step. The particular choice of the radius of analysis parameter 
is dictated to avoid any possible contamination by very off-axis residual particle events. 



^http : / / agile . asdc . asi . it 
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3.2. Results 

Figured] shows a Gaussian-smoothed intensity map (~ 10° x 08°) in Galactic coordinates 
integrated over the whole observing period, using the selections described in §3.11 The source 
detection significance is 19-a and the average 7-ray fiux above 100 MeV for the whole period 
is -Fe>ioomcV = (170 ± 13) x 10^^ photons cm~^ s~^, as derived from the AGILE Maximum 
Likelihood Code analysis. We note that the average 7-ray fiux computed over the three 
week campaign is lower than -pEMooMeV — (280 ± 40) x 10~® photons cm~^ s~^ , observed 
during the fiaring episode in July 2007, and computed during only a six day observation. 
Nevertheless, Figure [2] shows that the current average fiux is still higher than those observed 
during the EGRET era. The smaller errors on the AGILE November data with respect 
to the July data are due to both the higher statistics (323 versus 101 counts collected in 
November and in July, respectively), and to the more accurate calibration response files. 

Figure [3] shows the 7-ray light-curve at 1-day resolution for photons above 100 MeV. 
We note that 3C 454.3 is detected at a 3-0" level during almost the whole period on a 1-day 
timescale; this clearly indicates strong 7-ray fiaring activity. 

The average 7-ray fiux as well as the daily val ues of the 18 days were derived ac- 



cording to the 7-ray analysis procedure described in iMattox et al.l (119931 ). First, the en- 
tire period was analyzed to determine the diffuse gas parameters and then the source fiux 
density was estimated independently for each of the eighteen 1-day periods with the dif- 
fuse parameters fixed at the values obtained in the first step. Fitting the GRID fiuxes 
to a constant model (the weighted mean of the 1 -day average fiux v alues) yields F^td = 



1^186.3 ± 14.6) X 10 ^photons cm ^ s ^. Following iMcLaughlin et al.l (119961 ). we computed 



the V variability coefficient. The 2 -a upper limits (U L) were properly treated, assigning a 



value and a sigma equal to UL/2 (ITorres et al.ll200ll ). We obtain a = 36.7 for 17 de- 
grees of freedom (d.o.f.); therefore we can exclude that the fiuxes are constant at the 99.6% 
(~ 2.9cr) level, and we obtain a value for the variability coefficient V of 2.43. A value of 
V > 1 indicates that the source is variable within the observed period. 

Figure H] shows the average 7-ray spectrum derived over the entire observing period. 
The average spectrum was obtained by computing the 7-ray fiux in five energy bins over the 
entire observing period: 50 < E < 100 MeV, 100 < E < 200 MeV, 200 < E < 400 MeV, 
400 < E < 1000 MeV, and 1000 < E < 3000 MeV. We fit the data by means of a simple 
power-law model and restricted our fit to the most reliable energy range (100 MeV-1 GeV): 

/ -n^ \ -{1.73±0.16) 

F(E) = 3.61 X 10^^ I ph cm-2 s'^ MeV'^ (1) 



Unfortunately, the source was located substantially off-axis in the Super-AGILE FoV 
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during the whole observation period, resulting in a not particularly constraining upper limit 
flux, being as high as 1.13 x 10~^ photons cm~^ s~^ (50 mCrab). 



4. INTEGRAL data 

4.1. Data Reduction and Analysis 

The ESA INTEGRAL 7-ray Observatory, launched in 2002 October, carries three co- 
aligned code d mask telescopes. F or the purpose of this paper we refer to data from the IBIS 



1 



instrument (lUbertini et al.ll2003l ). sensitive in the energy range 15keV-10MeV and with a 



FoV of 29° X 29° , and in particular to the ISGRI lower energy detector layer. 

All the observations are organized into un-interrupted 2000-3600 s long science windows 
(SOW): light curves and spectra were extracted for each individual SCW. Wide-band spectra 
(from 17 to 150 keV) of the source were obtained using data from IBIS instrument. All the 
data were processed using the Off-line Scientific Analysis (OSA) version 7.0 software released 
by the INTEGRAL Scientific Data Centre. INTEGRAL data were analyzed using FTOOLS 
and XSPECll.3.2 in the Heasof t package (v.6.4). We assumed a single power law model to 
fit the IBIS data. 

Figure [5] shows the INTEGRAL/IBIS light curve in the energy range 20-50 keV ac- 
cumulated during the whole observation. The source does not show statistically significant 
flux variations. Figure [HI shows the INTEGRAL/IBIS spectra for revolution 623 (red circles), 
revolution 624 (blue triangles) and for the whole observation (black squares). 

Table [1] summarizes the INTEGRAL/IBIS spectral fit results. 



5. Swift data 



5.1. Data Reduction and Analysis 



The NASA Swift 'j-Taj Burst Mission (iGehrels et al.ll2004l ). launched in November 2004, 
has t hree co-aligned instruments: a coded- mask Burst Alert Telescop e (BAT; iBarthelmy et al 
15-150 keV), an X-ray Telesco pe fXRT: burrows et all boosl. 0.2-10 keV), and an Ul- 



2005 



traviolet/Optical Telescope (UVOT: iRoming et al 



2005 



170-600 nm) 



Swift data (Obs. ID 00031018) were collected by activating a Swift Cycle-3 Proposal 
(PI A. W. Chen) and by means of a dedicated ToO triggered by AGILE (PI S. Vercellone). 
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The XRT data were processed with standard procedures (xrtpipeline vO.11.6), adopting 
the standard filtering and screening criteria, and using FTOOLS in the Heasoft package 
(v.6.4). The source count rate was low during the whole campaign, thus we only considered 



photo n counting data (PC) and further selected XRT event grades 0-12 (IBurrows et al. 



20051 ). Swift/XKT data show an average count rate of > 0.5 counts s~^ and therefore pile-up 
correction was required. We extracted the source events from an annular source extraction 
region with an inner radius of 2-3 pixels (estimated case by case by means of the PSF fitting 
technique) and an outer radius of 30 pixels (1 pixel ~ 2'.'37). To account for the background, 
we also extracted events within a circular region centered on a region free from background 
sources and with radius of 80 pixels. Ancillary response files were generated with xrtmkarf , 
and account for different extraction regions, vignetting and PSF corrections. We used the 
spectral redistribution matrices vOlO in the Calibration Database maintained by HEASARC. 
Swift/XBT uncertainties are given at 90% confidence level for one interesting parameter (i.e., 
Ax^ = 2.71) unless otherwise stated. 

Figure [7] shows the 0.3-10 keV spectra for segment 001 to 006, where we summed seg- 
ments 003 and 004 in order to have similar statistics as the others: black circles (segm. 001), 
red squares (segm. 003+004), green upside-down triangles (segm. 005), and cyan stars 
(segm. 006). Segment 002 was not considered since only 1 s of data were recorded. We first fit 
Swift/XKT spectra with an absorbed power law model, named model A (wabs*zwabs (powerlaw) 
in XSPEC 11.3.2). Data were rebinned in order to have at least 20 counts per energy bin. 



The Galactic absorption was fixed to the value of N^''^ = 0.724 x 10^^ cm'^ flKalberla et al.l 



20051 ) . We considered, in addition to the Galactic absorption coeffi c ient, a n ext ra absorption 



compo nent, A^^5 following the results showed in iGhisellini et al.l (120071 ) and iRaiteri et al. 



(120071 ). A second spectral fit (model B) was performed considering a simple power law with 
the absorption component as a free parameter (wabs* (powerlaw)). Table [2] summarizes the 
most relevant spectral fit parameters. 

We note that the component in model A is not well constrained, while a simpler fit 
(model B) characterized by a single power law model with free absorption coefficient yields 
A^H values which are consistent within the uncertainties among dif ferent observations. We 
also check the possible presence of a double pow er law as reported inlRaiteri et al. J2008bh as 
follows. We fixed the Nu to the value derived bv lvillata etaP J2006h . (1.34±0.05) x 10^^ cm'^ 
(based on Chandra data), and fixed the hard power law index to the one we obtained by 
fitting the XRT data above 2keV. In most cases the harder power law component is not 
required, as its normalization is consistent with zero. 

Figure [8] shows the Swift/XRT photon index versus the 0.3-10 keV flux. Numbers 
beneath each point represent the observing segment. 
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We analysed Swift/BAT Survey data in order to study the hard X-ray emission of 
3C 454.3 and to investigate its evolution as a function of time. We selected two time windows, 
the first between 2005 A pril 01 and 2005 Sept ember 30 (when intense activity was recorded 
from the target, see e.g. iGiommi et al.ll2006l ). the second between 2007 June 01 and 2007 
December 31. We considered all BAT observations with 3C 454.3 in the FoV. After a careful 
data selectio n, based on backgro und rate, pointing stability and several other data quality 



criteria 
time of 



see 



Senziani et al.l 120071 ). we ended up with 4824 observations for a net exposure 
792 ks in 6 months in 2005 and ~ 624 ks in 7 months in 2007. 



Here we provide a brief overview of the BAT Survey data analysis procedure, while a 
detailed description will be addressed in a forthcoming paper (Belfiore et al., in preparation). 
Starting from Detector Plane Histogram (DPH) files. Detector Plane Images (DPI) in the 
20-60 keV and 60-100 keV energy ranges were generated and were cleaned from hot pixels 
and noisy detectors. Then, with the HEASOFT task hatfftimage each DPI was deconvolved 
to obtain a background- and source-subtracted sky image of the BAT FoV. For each sky 
image an appropriate effective exposure map (weighted on the coded fraction) was generated, 
accounting for possible Earth/Moon occultations. Then, such sky images were re-projected 
and stacked (weighting on effective exposure) to obtain monthly count rate maps, considering 
a small portion of the field around the target (a 3° x 3° region in local tangential projection 
coordinates, TAN). 

The count rates of 3C 454.3 were normalized to the Crab count rates and then converted 
to flux (ph cm~^ s~^) assuming for the Crab the canonical power law spectrum (photon index 
r = 2.15 and normalization of 10.4 ph cm~^ s~^ keV~^ at 1 keV) and for 3C 454.3 a power 
law spectrum with F = 1.7, averaging the instrument response over the field of view. The 
November 2007 net exposure time is ~ 106 ks for a flux in the 20-60 keV energy band of 
(1.07 ± 0.19) X 10-3 ph cm-2 s^^ 

Figure [9] shows the long-term Swift/BAT light curves in the 20-60 keV (bottom panel) 
and 60-100 keV (upper panel) energy range. The yellow vertical area marks the AGILE 
November campaign. The Swift/BAT flux is in good agreement with the flux derived from 



the whole INTEGRAL/IBIS campaign in the same energy range, -^^20*— 1 
ph cm" 



eOkeV 



1.02 X 10- 



-2 s-\ 



T he short-dashed line ma rks the epoch of the giant optical flare in 2005 



(IFuhrmann et al.ll2006l : IVillata et al.ll2006l ). when the hard X-ray flux was about twice higher 
than in November 2007. 
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6. Optical data 



6.1. WEBT Data Reduction and Analysis 

The Whole Earth Blazar Teles cope (WEBT) @ h as been monitoring 3C 454.3 si nce 



the exceptional 2004-2005 outburst fiVillata et al.l 120061 . 120071 : iRaiteri et al. 



20071. 12008al lb[). 

covering also the period of the AGILE observation. We refer to IRaiteri et al.l ( 2008al ) for 
a detailed presentation and discussion of the radio, mm, near-lR, optical and Swift/UVOT 
data collected, almost continously, during November 2007. 



6.2. REM Data Reduction and Analysis 



The photometric Optic al and Near Infra red (NIR) observations were carried out with 
the Rapid Eye Mount (REM. IZerbi et al.ll2004l ). a robotic telescope located at the ESO Cerro 
La Silla observatory (Chile). The REM telescope has a Ritchey-Chretien configuration with 
a 60 cm f/2.2 primary and an overall f/8 focal ratio in a fast moving alt-azimuth mount pro- 
viding two stable Nasmyth focal stations. At one of the two foci th e telescope simultane ously 
feeds, by means of a dich roic, two cameras: REMIR for the NIR (jConconi et al.l 120041 ). and 
ROSS (ITosti et al.ll2004l ) for the optical. Both the cameras have a field of view of 10' x 10' 
and imaging capabilities with the usual NIR {z\ J, H and K) and Johnson-Cousins V , R, 
and / filters. All raw optical/NIR frames obtained with REM Telescopes, were corrected 
for dark, bias and fiat field. Instrumental magnitudes were obtained via aperture photom- 
etry using GAIA0. Calibration of the optical source magnitude was o btained by differentia l 
pho tometry with r espec t to the comparison stars sequence reported by lFiorucci et al.l (119981 ) 
Raiteri et al. Nl998h. For th e NIR calibration we used the comparison sequence reported 



and 



by I Gonzalez- Perez et al.l (120011 ). Both REMIR and ROSS instruments were used in order 
to obtain nearly simultaneous data and to study the spectral behavior 3C 454.3 at different 
levels of flux. 



7. Simultaneous data analysis 

Figure [TU] shows the simultaneous light curves acquired during the period 2007 Novem- 
ber 6-December 3. Black circles represent AGILE/GRID data (30 MeV-50 GeV); red tri- 



■^http : //www . oato . inaf . it/blazar/webt, see e.g. IVillata et al.l (j2004l ). 
'*littp : / / star-www . dur . ac . uk/ ^pdraper/gaia/gaia . html 
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angles represent INTEGRAL/IBIS data (20-200 keV); blue pentagons represent Swift/X.KT 
data (0.3-10 keV); c yan-solid and green-open squares represent i?-band WEBT and REM 
(iRaiteri et al.ll2008al ) data, respectively. The yellow areas mark the periods PI and P2 dur- 
ing which we compute the simultaneous spectral energy distributions, and corresponding to 
higher 7-ray flux levels. We note that during the period PI the optical flux shows intense 
variability, reaching a relative maximum on the last day of the 7-ray day-by-day sampling. 
Moreover, an optical flare as intense as the one on MJD~54420 occurred at the end of the 
AGILE observations (MJD~54435.5). The three 7-ray data points show no particular flaring 
activity, though the daily flux remained quite high, ~ 200 x lO^^photons cm~^ s~^. A de- 
tailed discussion on the correlation between the optical and 7-ray data during December 2007 
wiU be presented in the forthcoming Paper II (Donnarumma et al. 2008, in preparation). 

We investigated th e expected 7-optical flux correlation by means of the discrete correla - 
tion function (DCF; see lEdelson fc KroliklllQSSi : iHufnagel fc Bregmanlll992l : |Petersonll200ll ). 
The result is shown in Figure [TTl The DCF peak occurred at r = 0, and its value is ~ 0.5. 
This indicates a moderate correlation, with no signiflcant time delay between the 7-ray and 
optical flux variations. A minor peak at r = —5 days comes from establishing a connection 
between the optical flare at MJD ~ 54420 and the 7 high flux at MJD ~ 54425. 

Figure [12] shows the spectral energy distribution for the period PI, MJD 54417.5-54420.5 
(see Figure [To]). Filled squares represent the AGILE/GRID data in the energy range 100- 
1000 MeV; small flUed circles represent Swift/XKT data in the ene rgy range 0.3-10 keV 
(segment 001); open symbols represent radio to UV data taken from iRaiteri et al.l (l2008a] ). 
corresponding to MJD 54420, when all the WEBT UBVRI bands were available, as well as 
Swift/VVOT data. 

Figure [T^ shows the spectral energy distribution for the period P2, MJD 54423.5-54426.5 
(see Figure [TUl) . Filled squares represent the AGILE/GRID data in the energy range 100- 
1000 MeV; flUed triangles represent INTEGRAL/IBIS data in the energy range 17-150 keV 
(orbits 623-1-624); small flUed circles represent Swift/XBT data in the energy range 0.3- 
10 keV (segrnents 00 3, 004, and 005); open symbols represent radio to UV data taken from 
Raiteri et al.l (l2008al ). corresponding to MJD 54425. 



A brief discussion of the modeling of both SEDs is presented in the Section [8] 



8. Discussion 



The long-term 7-ray activity of 3C 454.3 is one of most interesting discoveries achieved 
by AGILE during its flrst 6 months of observations. The source was already detected in high 
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state in July 2007 during a 1-week AGILE ToO triggered by an intense optical flare detected 
by the WEBT. During that period, the source reached its highest intensity level, with an 
average 7-ray flux of -FE>iooMeV = (280 ±40) x 10~^ photons cm~^ s~^. In November 2007, 
3C 454.3 showed prominent and prolonged 7-ray activity, with flaring episodes on a timescale 
of a few days and an average 7-ray flux of Fe> looMeV = (170 ± 13) x 10~^ photons cm~^ s~^. 
This renewed activity triggered observations at different frequencies, allowing us to obtain 
an almost simultaneous SED coverage on 14 decades in energies. We dubbed 3C 454.3 as 
crazy diamond because of its dramatic variability at high energies revealed during the flrst 
half of the AGILE Observing Cycle-1. It has become clear that this source is playing the 
same role for AGILE as 3C 279 had for EGRET. The 3 C 454.3 stron g variab ility has also 



a clear signature at lower frequencies. As reported in iRaiteri et al.l (l2008al ). during the 
AGILE observation, on MJD 54425 the source showed an extremely variable behavior in 
the R band, with a brightening of 0.33 mag in 2.3 hours. In the same paper, the authors 
report other episodes of fast variability with flux variations of several tenths of mag in a 
few hours. Moreover, while in July 2007 3C 454.3 exhibited its most intense optical flare 
but with a very moderate degree of 7-ray flux variability on a day-by-day time-scale, during 
the November 2007 campaign (see Figure [3]), we note a signiflcant 7-ray flux variability on 
short time-scales with at least two distinct flaring episodes (PI and P2). It is interesting 
to note that in the optical band 3C 454.3 also seems to display more rapid flares during the 
fall-winter 2007 WE BT campaign than those that occurred during the July 2007 monitoring 



flRaiteri et al.ll2008bf ). 



We compared the spectral properties of higher-state periods, PI and P2, with two 
lower-state periods, PJowl and P_low2, chosen of the same duration as PI and P2 and 
corresponding to MJD 54414.5-54417.7 and MJD 54420.5-54423.5, respectively. Figure [14 
shows the AGILE/GRID spectra for periods PI (red square), P2 (green star), PlJowl (black 
circle), P2Jow2 (blue upside triangle). The July 2007 spectrum is also shown (cyan upside 
down triangle). Although the statistics accumulated in only 4 days does not allow us to 
obtain a robust flt of the data. Figure [H] shows no clear spectral differences among different 
source intensity levels. 

The correlation between the flux level and the sp ectral slope was ext e nsive ly studied by 



means of the analysis of the EGRET data. Recently, iNandikotkur et al.l (120071 ) have shown 
that there is no homogeneous behavior among EGRET blazars. Although they consider 
long-term spectral dependence on flux rather than short-t erm as in our case, our fl ndings 



are in agreement with their results on 3C 454.3. Figure 3 in INandikotkur et al.l (l2007l ) shows 



no spectral variation despite a flux variation of about a factor of four. 

Different emission mechanisms can be invoked to explain the 7-ray emission and the 
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different spectral states. In the leptonic scenario, the low-frequency peak is interpreted as 
synchrotron radiation from high-energy electrons in the relativistic jet, while the high-energy 
peak can be produced b y IC on differerit flav o rs of seed photons. In th e synchrotron self- 
Compton [SSC] model (iGhisellini et al.l Il985l . iBloom &: Marscherl Il996l ) the seed photons 
come from the jet itself. Alternatively, the seed photons can b e those of the accre tion 
disk [external Compton scattering of direct disk radiation, ECD, iDermer et al.l (119921 )] or 
those of the broad- line region (BLR) clouds [external Compton scattering from clouds, ECC, 
Sikora et al.l (Il994j )]. The target seed photons can also be those produced by the infrared (IR) 
dust torus surroun ding the nucleus [external Compton scattering from IR dust, ERC(IR), 
Sikora et all J2002h ]. 



Blazejowski et al.ll2000l showed that the ERC(IR) emission peak is at lower frequencies 



(soft 7-ray ), more suitable for MeV blazars, while our SEDs show a peak for the EC com- 
ponent in the GeV region of the spectrum. The average photon index during this November 
2007 campaign (r A OTT.R = 1-73 ± 0.16) is harder than the time-averaged one reported in 
Nandikotkur et al.l (120071 ) (Eegret = 2.22 ± 0.06) for EGRET. During intense 7-ray flares, 
the ECC and ECD processes play a major role and the softness or the har dness of the re- 



sultin g spectrum is controlled by the dominant component, as illustrated in iHartman et al. 



( 120011 ) for 3C 279. In the case of 3C 454.3, the ECC component seems to play a major role, 
as we will show from the SED modeling. 

We fit the SEDs for the PI and P2 gamma-ray flaring episodes by means of a one- zone 
leptonic model, considering the contributions from SSC and from external seed photons 
originating both from the accretion disk and from the BLR. The emission along the jet is 
assumed to be produced in a spherical blob with comoving radius R by accelerated electrons 
characterized by a comoving broken power law energy density distribution of the form. 



Ml) 



(2) 



(7/7b)°i + (7/7b)"^ ' 

where 7 is the electron Lorentz factor assumed to vary between 10<7<1.5xlO^, ai and 
Oh are the pre- and post-break electron distribution spectral indices, respectively, and 7b 
is the break energy Lorentz factor. We assume that the blob contains a random average 
magnetic field B and that it moves with a bulk Lorentz Factor F at an angle Gq with respect 
to the line of sight. The relativistic Doppler factor is then 6 = [F (1 — /? cosGq)]"^, where /3 
is the usual blob bulk speed in units of the speed of light. 

Our modelling of the 3C 454.3 high-energy emission is based on an Inverse Compton 
model with two main sources of external target photons: (1) an accretion disk characterized 
by a blackbody spectrum peaking in the UV with a bolometric luminosity for a IC- 
scattering blob at a distance L from the central part of the disk; (2) a Broad Line Region 
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with a spec trum peaking in the V band and assumed to reproces s 10% of the irradiating 
continuum JXavecchio fc Ghiselhnilboosl : iRaiteri et aLlboOTI . boOSbh . 



These two regions contribute to the ECD and the ECC, respectively, and it is interesting 
to test the relative importance of the two components that can be emitted by the relativistic 
jet of 3C 454.3 under different conditions. A complete theoretical analysis of the model, and 
of the interplay among the different parameters is beyond the scope of this paper and it will 
be presented in a forthcoming paper. We summarize here the main results of our best model 
characterized by an interesting set of physical parameters. 

Table [3] shows the best-fit parameters of our modelling of the flaring state SEDs corre- 
sponding to the PI and P2 phase of Figure [TOl Our best fit parameters values are: 5 ~ 10 G, 
r = 8.4 , 00 = 2.6° , and r = 0.05 pc where r is the distance between the accretion disk and 
the emitting region. In both Figures [12] and [121 the dotted, dashed, and dot-dashed lines 
represent the contributions of the accretion disk blackbody, the external Compton on the 
disk radiation and the external Compton on the broad line region radiation, respectively. 

We note that during both the PI and P2 episodes, the ECD contribution can account 
for the soft and hard X-ray portion of the spectrum, which show a moderate, if any, time 
variability. However, we note that the ECD component alone cannot account for the hardness 
of the 7-ray spectrum. We therefore argue that in the AGILE energy band a dominant 
contribution from ECC seems to provide a better fit of the data during the gamma-ray 
flaring states PI and P2. Moreover, Table [3] shows that the data of both SEDs can be fit by 
very similar model parameters. We note, however, that the high energy part of the electron 
energy distribution appears to be softer dur i ng th e P2 episode as compared to the electron 
distribution of the PI flare. iHartman et al.l (12001! ) find, for 3C 279, a relevant contribution 
of the SSC component in the X-rays - soft 7-ray bands. However, their average values for 
7i and 72 are much higher (a factor of 5-10) than ours, resulting in an increase of the ratio 
SSC/Sync oc 7^ . Thus, for the Hartman et al. choice of parameters, the SSC contribution 
becomes relevant at higher frequencies and of the same order of the ECD contribution. 



Our results can be compared with those obtained by jjorstad et al.l (120051 ). By means of 
total and polarized images obtained at the Very Long Baseline Array (VLBA) at 7mm, they 
were able to compute the global parameters of the source jet, estimating (P) = 15.6 ± 2.2 , 
{S) = 24.6 ± 4.5 , ^ = (0.8 ± 0.2)° , and (Go) = (1.3 ± 1.2)° , where (P) and {S) are the 
average Lorentz and Doppler factors, respectively, while 6 and (Gq) are the intrinsic half 
opening angle of the jet and the angle between the jet axis and the line of sight, respectively. 
We note, however, that the jet parameters derived so far were obtained by means of data 
collected in earlier epochs with respect to our observations and refer to average values of 
different jet components. 
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The energetics of 3C 454.3 can be computed by estimating the isotropic luminosity in 
the 7-ray band, L]^°, and comparing it with the Eddington, the bolometric, and the particle 
injection luminosities. For a given source with redshift z, the isotropic emitted luminosity 
in the energy band e is defined as 

i(,). = _i£M£)_, (3) 

where, in our case, e is the 7-ray energy band with -Emin = 100 MeV and -Emax = 10 GeV, 
a is the 7-ray energy spectral index, F(z/) oc z/~" is the energy differential fiux, F = 
rEnia.^/h pt^\^^ -g ^YiQ flux in the 7-ray band, and the luminosity distance is given by 

d,{zi,Z2) = {I + Z2)' X r[E{z)]-'dz, (4) 

J- + ^2 Jzi 

where zi = 0, Z2 = Zsrc and 

E{z) = ^ynuii + z)^ + (1 - - ^^a)(1 + z)^ + , (5) 

where Hq is the Hubble constant, Qm and are the contribution of the matter and of 
the cosmological constant, respectively, to the density parameter. Hereafter, we assume 
Hq = 70kms~^Mpc~^, Qm = 0.3, and ^a = 0.7. Using the observed average 7-ray fiux, we 
obtain L]^" = 3.9 x 10^^ erg s-\ 

Moreover, from the values quoted in Table [3] and from Equation [2] we can compute the 
particle injection luminosity, Ljnj , obtaining : 

Linj =7iR^T^c j [c/7 me 0^7 n(7)] = 3 x 10^^ erg s"^ . (6) 



Assuming for 3C 454.3 a black hole mass Mbh = 4.4 x lO^M© fiGu et all 120011 ) . we 



obtain an Eddington luminosity of the order of L-^aa = 5.7 x lO^^ erg s ^ to be comp ared 



with the bolometric luminosity Lboi = 1-9 x lO^^erg s ^ reported in IWoo &: Urryl (120021 ). 



We obtain, therefore, that the source energetic is comparable to the value obtained by 



Tavecchio et al.l (120071 ) for the power of the inner portion of jet, few x 10^^ erg s ^. 



9. Conclusions 



The AGILE mission is particularly suited to monitor a large number of potential 7-ray 
sources. The AGILE pointings during the month of November 2007, despite being centered 



-17- 



approximately in the Cygnus region of the Galactic plane, revealed the very prominent 7-ray 
activity of the blazar 3C 454.3. The AGILE detection of this blazar prompted a series of 
important multiwavelength observations. The electromagnetic emission of 3C 454.3 could 
be determined with unprecedented coverage over 14 orders of magnitude in energy during a 
period that included a substantial fraction of the months of November and December, 2007. 
Results of the AGILE data and related multifrequency campaign carried out in December, 
2007, will be presented in a forthcoming paper. 

We reported in this paper the main results of our AGILE and associated multifrequency 
campaigns during the month of November, 2007. Our results can be summarized as follows: 

1. The 7-ray emission from 3C 454.3 dominated the whole extragalactic sky as monitored 
by AGILE during its first year of scientific operations. 

2. Our 7-ray data show remarkable variability on a daily timescale for a Flat Spectrum 
Radio Quasar. 

3. Emission in the optical range appears to be correlated with that at 7-ray energies 
above 100 MeV. 

4. Variability in the soft and hard X-ray range is less sensitive to the short-timescale 
variations of the optical flux. 

5. The average 7-ray spectrum during the whole campaign is substantially harder than 
that reported in previous observations. 

6. We determined the SEDs for episodes of relatively high 7-ray emission. 

7. Our results support the idea that the dominant emission mechanism in 7-ray energy 
band is the inverse Compton scattering of external photons from the BLR clouds 
scattering off the relativistic electrons in the jet. 

We thank the Referee for his/her very prompt and constructive comments. The AGILE 
Mission is funded by the Italian Space Agency (ASI) with scientific and programmatic par- 
ticipation by the Italian Institute of Astrophysics (INAF) and the Italian Institute of Nuclear 
Physics (INFN). This investigation was carried out with partial support under ASI contract 
N. 1/089/06/1. This work is partly based on data taken and assembled by the WEBT col- 
laboration and stored in the WEBT archive at the Osservatorio Astronomico di Torino - 
INAF (http://www.to.astro.it/blazars/webt/). We thank the Swift and INTEGRAL 
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planners. MTF, AB and PU acknowledge ASI/INAF contract N. 023/05/0. 
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Table 1. INTEGRAL/IBIS spectral fit results. 



Rev. 


r 


(d-o.f.) 


Flux ^ 


623 

624 

623+624 


1 70+0.33 

-■-•'^-0.30 
1 71+0.41 
1- ' 1-0.36 
1 7C+0.25 
-I- '^-0.23 


1.21 (11) 

0.54 (11) 
0.89 (11) 


1.52 

1.42 
1.49 


(")Flux 


in the 20- 


200 keV band 


in units 



of 10 ^photons cm ^ s ^obtained from the 
spectral fits. 



Table 2. Swift/XRT spectral fit results. 



Obs.^ 


(10^2 cm2) 


r 


Xljid.ol.] 


1 Fo.3-10 (obs.) 
photons cm~2 s~^ 


-^2-10 
photons 


(obs.) 

— 2 —1 

cm s 





Model A: single power law with Galactic+intrinsic absorption 



001 0.0724+(0.06+°:°^) 1.56;°;°^ 1.173/(108) 1.45 x lO'^ 0.52 x 10"^ 

003+004 0.0724+(0.20lg:iJ) 1.681^;^^ 0.933/(53) 1.41 x 10"^ 0.51 x 10"^ 

005 0.0724+(0.23l[5:J^) l.GOtoH 0.900/(46) 1.43 x 10"^ 0.56 x 10"^ 

006 0.0724+(0.25i^:?^) 1.681^;^^ 1.089/(80) 1.28 x 10"^ 0.48 x 10"^ 



Model B: single power law with free absorption 



001 0.10t[J:[!^ 1-561^7 1.100/(131) 1.39 x 10"^ 0.51 x lO'^ 

003+004 0.15l[J:|J^ l-73^o;i2 0.921/(53) 1.40 x 10"^ 0.50 x lO'^ 

005 0-l5-om l-66+o:n 1.070/(58) 1.35 x lO'^ 0.52 x 10"^ 

006 0.16^°:°^ 1.124/(94) 1.23 x 10-^ 0.46 x lO'^ 



'^"•'ast three digits of observation number 00031018 
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Table 3. Input parameters for the model of PI and P2 SEDs. See text for details. 



Parameter 


SED PI 


SED P2 


Units 


ai 


2.1 


2.2 




ah 


4.5 


5.0 




Tmin 


10 


10 




7b 


500 


500 




K 


14 


12 


cm~^ 


R 


35 


35 


10^^ cm 


B 


10 


8 


G 


5 


14.64 


14.64 




La 


5 


5 


10^6 erg s"^ 


r 


0.05 


0.05 


pc 


©0 


2.6 


2.6 


degrees 


r 


8.4 


8.4 





Fig. 1. — Gaussian-smoothed intensity map (~ 10° x 08°) in Galactic coordinates integrated 
over the whole observing period (2007 November 10 12:17 UT - 2007 December 01 11:39 
UT). The cross symbol is located at the 3C 454.3 radio coordinates. [See the electronic 
edition of the Journal for a color version of this figure] 
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Fig. 2. — EGRET (triangles) and AGILE-GRID (square ar id circle) gamrn a -ray l ight curve 
in units of 10~^ photons cm~^ s~^. EGRET data are from iHartman et al.l (jl999l ). AGILE 
July 2007 data are from V08. 
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Fig. 3. — AGILE-GRID 7-ray light-curve at ~ 1-day resolution for E>100 MeV in units of 
10~^ photons cm~^ s~^. The downward arrows represent 2-a upper-limits. The dashed line 
represents the weighted mean flux. 
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Fig. 4. — AGILE-GRID average 7-ray spectrum. Three energy bins were considered: 100 < 
E < 200 MeV, 200 < E < 400 MeV, 400 < E < 1000 MeV. The red dashed line represents 
the best-fit power law model. [See the electronic edition of the Journal for a color version 
of this figure] 
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Fig. 5. — INTEGRAL/IBIS light curve in the energy range 20-50 keV accumulated during 
the whole observation. The dashed line represents a fit with a constant model. [See the 
electronic edition of the Journal for a color version of this figure] 
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Fig. 6. — INTEGRAL/IBIS spectra for revolution 623 (red circles), revolution 624 (blue 
triangles) and for the whole observation (black squares). [See the electronic edition of the 
Journal for a color version of this figure] 
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Fig. 7. — Swift/XBT 0.3-10 keV spectra for segment 001 to 006: black circles (segm. 001), 
red squares (segm. 003+004), green upside-down triangles (segm. 005), and cyan stars 
(segm. 006). [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 8. — Swift/XRT photon index versus the 0.3-10 keV flux. Numbers beneath each point 
represent the observing segment. Segment number 2 is missing since only 1 sec of data were 
recorded. 
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Fig. 9. — Long-term Swift/BAT light curves in the 20-60 keV (bottom panel) and 60- 
100 keV (upper panel) energy range. The yellow vertical area marks the AGILE November 
campaign. The short-dashed line marks the epoch of the giant optical flare in 2005. [See the 
electronic edition of the Journal for a color version of this figure.] 
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Fig. 10. — Simultaneous light curves acquired during the period 2007 November 6-December 
3. Black circles represent AGILE/GRID data (30 MeV-50 GeV); red triangles represent IN- 
TEGRAL/IBIS data (20-200 keV); blue pentagons represent Swift/XR T data fO.3-10 keV) : 
cyan-solid and green-open squares represent /2-band WEBT and REM ( iRaiteri et al.ll2008al ) 
data, respectively. The yellow areas mark the periods PI and P2 during which we compute 
the simultaneous spectral energy distributions. [See the electronic edition of the Journal for 
a color version of this figure.] 
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Fig. 11. — Discrete correlation function between the 7-ray and optical fluxes. The optical 
data have previously been binned over 12 hours to smooth the intranight variations. The 
DCF peak suggests a mild correlation with no time delay. [See the electronic edition of the 
Journal for a color version of this figure] 
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Fig. 12. — Spectral energy distribution for the period PI, MJD 54417.5-54420.5 (see Fig- 
ure [10]). Filled squares represent the AGILE/GRID data in the energy range 100-1000 MeV; 
small filled circles represent Swift/XBT data in the en ergy range 0.3-10 keV (segment 001); 
open symbols represent radio to UV data taken from iRaiteri et al.l (j2008al ). corresponding 
to MJD 54420. The dotted, dashed, dot-dashed, and triple-dot dashed hues represent the 
accretion disk, the external Compton on the disk radiation, the external Compton on broad 
line region radiation, and the SSC contributions, respectively. 
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Fig. 13. — Spectral energy distribution for the period P2, MJD 54423.5-54426.5 (see Fig- 
ure [TO]). Filled squares represent the AGILE/GRID data in the energy range 100-1000 MeV; 
filled triangles represent INTEGRAL/IBIS data in the energy range 20-200 keV (orbits 
623+624); small filled circles represent Swift/X.KT data in the energy range 0.3 -10 keV (seg- 



ments 03, 004, and 005); open symbols represent radio to UV data taken from lRaiteri et al. 



(l2008al ). corresponding to MJD 54425. The dotted, dashed, dot-dashed, and the triple-dot 
dashed lines represent the accretion disk, the external Compton on the disk radiation, the 
external Compton on broad line region radiation, and the SSC contributions, respectively. 
[See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 14. — AGILE/GRID spectra for periods PI (red squares), P2 (green stars), PlJowl 
(black circles), P2_low2 (blue upside triangles). The July 2007 spectrum is also shown (cyan 
upside down triangles). 



